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1. INTRODUCTION 


The planning of the Firsf GARP Globa! Experiment includes a subsystem 
of superpressure balloons (WMO-ICSU, 1969, p. 10). Earth orbiting satellites v/ill 
track the motions of these balloons, and communicate with sensors on them, and 
thereby provide observations of the wind, temperature, pressure and geopotential of 
the air at the position of each balloon (cf. Lally and Lichfield, 1969; Morel, 1969). 

A preliminary test of a satellite-balloon system of this kind, on a Hemi- 
spheric-scale, Is the "Hole Experiment", sponsored by the French Centre National 
des Etudes Spatiales in cooperation with the Planetary Radiations Branch of the Lab- 
oratory for Atmospheric and Biological Sciences of the U.S. National Aeronautics 
and Space Administration. 

The Eole Experiment Is scheduled for the Southern Hemisphere summer of 
1970-71 . According to the present plans, after the tracking and interrogation satellite 
has been placed in its polar orbit, 10 balloons will be released per day at each of the 
following launching sites in Argentina: 

A. Logo Fagnano .54.5°S 67.5°W, 

B. Neuquen 38.-6*^S 68.1°W, 

C. Mendoza 32.5°S 68.5°W, 

until a lota! of 500 balloons are in the air. These balloons will float on the surface 
of constant air density, which in the Standard Atmosphere has the pressure of 200 mb. 
The balloons are expected to have a 90 day mean lifetime. As individual balloons 
fail, the total number of 500 balloons will be maintained in the air by additional 
launchings, until all of the 750 balloons provided for the experiment have been used. 
Any balloon will be automatically cut down if it crosses into the Northern Hemisphere. 



2 . 


In the study presented in this report, we have attempted to anticipate 
the performance of the Eole system by simulating the trajectories of the balloons 
for a 130 day period. For this prupose we use the winds produced in a numerical 
simulation of the genera! circulation of the global atmosphere. 


2. NUMERICAL SIMULATION OF THE GLOBAL WIND FIELD 


The winds used in this study were obtained by numerical integration 
of the primitive equations with a two“levei version of the general circulation model 
described by Arakawa, Mintz and Katayama (1 970). In this mode! , the atmosphere 
between the earth's surface and 200 mb pressure is divided into two layers of equal 
mass. The time-dependent variables are the horizontal velocities and temperatures 
in the middle of the two layers, the surface pressure, the water vapor in the lower 
layer, and the water stored in the ground. The earth's surface is specified as open 
ocean, ice covered ocean, and bare land or snow covered land. The height of the 
land follows the large scale mountain systems. The temperature of the ocean surface 
and the distribution of the sea ice are prescribed as their observed annual mean values. 

The grid points are at every 5*^ of longitude and 4° of latitude, from 
pole to pole. The finite space-difference scheme maintains the integral constraints 
of physical importance, and the scheme suppresses any computational noise in the 
pressure and temperature fields which is not under control of the geostrophic adjust- 
ment process. The time step is 6 minutes, using a combination of the Heun and the 
Matsuno time-difference schemes. 
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Kinetic energy dissipation in the model is only by surface friction 
and friction due to the vertical eddy stress betv/een the two wind levels. The 
heating by absorption of solar radiation depends on the calculated water vapor 
and diagnostically determined clouds. The long wave radiative cooling depends 
on the water vapor and clouds and on the temperatures of the air and earth's surface. 
The sensible heat transfer, at the air-ground interface, depends on the temperature 
difference between the air and the ground. Over bare land and over snow and ice 
surfaces, the ground temperature is calculated from the surface energy budget. 

Latent heat release is obtained from the calculated precipation. 

There are two kinds of precipitation in the model . One is large scale, 
which occurs when the lapse rate is absolutely stable and the grid scale humidity 
exceeds its saturation value. The other kind is cumulus-convective precipitation. 
The cumulus-convection is of three forms; low-level, middle-level and penetrative 
convection, depending on the humidity and the vertical distribution of temperature. 
The low-level convection only transfers sensible heat and water vapor from the 
ground to the air. It forms low level convective clouds which influence the short 
and long-wave radiative transfers, but all the water in these clouds evaporates into 
the air. The middle-level and peneirative convection not only transfer heat and 
water vapor, but also release latent heat of condensation. 

Evaporation from the surface depends on the wind speed and on the 
difference between the vapor pressure of the air and the vapor pressure of the surface 
For open ocean, ice and snow, the surface vapor pressure is its saturation value. 

For bare land, the surface vapor pressure depends on the wetness of the ground. The 
ground wetness is calculated from the precipitation, evaporation and runoff. 

he initial state in this numerical simulation of the global atmospheric 
circulation, was a state of zero winds (atmosphere and earth in solid rotation); 
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fhree-dimensional isothermalcy^ at 250°K;a constant sea level pressure, of one 
atmosphere; no water~vapor in the air; and dry ground. 

The v/ind and mass fields used to calculate the balloon trajectories in 
this numerical simulation of the Eoie Experiment, are from the 130 day period, day 
341 to day 471 of the numerical integration, during which time the sun changed its 
declination and distance from that of 1 December to 1 1 April. 


3. CALCULATION OF BALLOON TRAJECTORIES 


In the numerical simulation; we calculate the trajectories of the balloons 
which float on the surface of constant air density, ~ 0.329 kg/m , which is the 
air density at 200 mb in the Standard Atmosphere. 


From the two levels at which the winds and temperatures are computed in 
the model, the wind components are extrapolated linearly with pressure to the density 
surface, p^, where they are used to transport the balloons. 

Although the lime-dependent fields In the numerical simulation of the 
global circulation were calculated- in 6 minute steps, for practical reasons these fields 
were recorded on magnetic tapes only at 6 hour Intei-vals. For the balloon trajectory 
calculations the 6 hourly winds were linearly interpolated to 1 .2 hour intervals. 
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Horizonfxj! Inteqaolat-ion, from riie grid points at which the winds were 
computed (every 5° of longitude and 4° of latitude) to the locations required for the 
trafectory calculation, was by linear Interpolation in longitude and latitude. 

After the winds were Interpolated vertically, in time, and horizontcrlly, 
to the locations of every balloon, 1 ,2 hour forward displacements were taken along 
the great circle paths given by the inteipolated winds. At the endpoint positions of 
the 1 .2 hour displacements, the time and space interpolations of the wind were 
repeated, and new 1 .2 hour forward displacements calculated. The end points of 
these second displacements were then connected to the starting points of the initial 
displacements by great circle segments, and the median points of these segments were 
taken as the positions of the balloons at the end of the first 1 .2 hour time step. This 
space- and time-centered procedure for computing tra|ectories on a spherical surface 
corresponds to the method given by Henrici (1962) for calculations on a plane, with 
his equation 3-47 and the constant, a = 0.5. 

Iteration of the above procedure in 1 .2 hourly steps, for the 500 balloons, 
and for 130 days, produced the distribution of balloons shown in the next section. 

In the numerical simulation, the balloons were introduced directly on 
the density surface, p^ , at intervals of 2.4 hours, at the location of each of the three 
launching stations, until 500 balloons were in the air. The first three balloons were 
introduced at time 0.0 GMT, December 1 , year 1 , In the mode! atmosphere; and the 
last two balloons at 14.4 GMT, December 17, year 1 . The balloons are identified 
by the lettei^ A,B,C according to their launch sites: A = Logo Fagnano; B=Neuquen; 
and C=Mendoza. 
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4. HEMISPHERIC MAPS OF THE BALLOON DISTRIBUriONS 


The calculated distributions of the balloons are shown/ in figures 
through 28, on polar equal area maps of the Southern Hemisphere. Because the 
positioning of the balloons is constrained in the figures by the resolution of the 
line-printer/ wherever two balloons fall within the same printing location the 
second balloon (and, more rarely / others too) is printed in the nearest vacant 
printing space at one or the other side/ or belov// the true printing location. 

On day 0 (hour 0.0), there are only three balloons shown on the maps,' 
one at each of the three launch sites. 

At day 2.5 there are 78 balloons in the air, and they show westerly 
steakline waves, of variable amplitudes. 

At day 5 there are 153 balloons, and the streaklines are still relatively 
well defined by balloons B and C. But balloons A, those launched at Logo Fagnano, 
have mostly separated into two loose clusters. 

By day 10 balloons B and C form a ring around the globe, with about 
7 small amplitude v/aves in the streaklines. Some of the balloons A are now In the 
polar region. By this time, 303 balloons are in the air. 

On day 15 there are 453 balloons in the air. But there are now fewer 
balloons in the polar region. 

By day 20, 500 balloons are in the air. These same balloons remain in 
the air for the remainder of the 130 days, as we do not simulate any balloon failures 
and replacements. 



UO^VCNIC Smr^CE KC/M««). ^ITUillQH iFTEA 0«0 OAVS 



Fig. 1. 


Balloon positions, on a polar equal area projection, at day 0 
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tsorvCHir su«f«fE no.o.ijo siTu»Tt(it< *ftf» 2.5 ni»5 



Fig. 2. Balloon positions, on a polar equal area proiection, at day 2.5. 


ISOPyC^tC SUdFACE ilTUftTlON 6FTEA ^,9 OAlfS 



Fig. 3. 


Balloon positions, on a polar equal area profection, at day 5. 
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iSOPYCNIC KO*0«S29 KG/NOaS* SITUAtIDK APTIft 10*0 9AYS 



Fig. 4. Balloon positions/ on a polar equal area projection/ at day 10. 


!$DPYCN1C SUFFICE R0«6*32^ SITUATION AFTER 


15*0 OAVS 



Fig. 5 


Bailoon positions, on a polar equal area projection, at day 15 


ISQ^YCNSC SUftFACC R0-0*52«» KG/»1«*3* SITUATlOK 4FT60 


20»0 C4YS 
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Balloon positions/ on a polar equal area projection/ at day 20 



ISOPTCmiii SITUATIttH *^T6« 2#*0 OtVS 
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Balloon positions/ on o polar equal area projection/ at day 



isnPYCHic sysFAce rD’O.iJ'J hc/h**3, 


StTUiTIQH AFFER 


30.0 DAYS 



Fig. 8 


Balloon positions, on a polar equal area protection, at day 



SITUATtOK AfTER 
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UDPrCNIC SUAFACP R0-0»52«* 


39.0 DAYS 



Fig. 9 


Balloon posi^^ons;. on a polar equal area projecHon, at day 35 
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TSO^VCHIC SURFACE UO*O.J2'» KC/««*3. SITUATION AFTER AO.O ’JAYS 



Fig. 10. Bolloon positions, on a poiar equal- area projection, at day 40 




ISOPYCNIC SUftfACe ao-o.j?*? KG/«**5. SITUATinn aftcr 
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45.0 DAYS 



Fig. 11 


Balloon positions, on a polar equation area projection/ at day 
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]SO?TCHtC SUitFiCE K0>0«329 K0/n««3* SITUATIDH jfTEA 50»0 OiYS 



Fig. 12. Balloon posiHons;. on a polar equal area projecHon, at day 50. 
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tSO^TCNIC SU«F4CE K0/«**5. SITUJITIOH AFTER 55.0 DAYS 



Fig. 13. Balloon posiHons, on a polar equal area projection, at day 55 


ISnPYCNIC surface R0»0.3^« SlUJATtON AFTER 


60^0 DAYS 



ig. 14. Balloon positions, on a polar equal area projection, at day 
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ISOfYCNlC StlKPACC Rn*0,A2<I SlTU'tTKm iFTFR 6S.0 DAYS 



Fig. 15. Balloon positions, on a polar equal area projection, at day 65 
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l$0?VCNtC SURFACE R0«0,5?^ KC/K**3« S;TUAT!0« AFT£R 70.0 OiYS 



Fig. 16. Balloon positions;^ on a polar equal'area projection, at day 70 
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iso^ycwir su‘iFAre pn»o>?^'9 kc/m**?. situation after t*j,o oj.ys 



Fig. 17 Balloon positions/ on a polar equal area projection/ at day 75 
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* tSOPYCKlC SURFACE Rn«0.329 SITUATION AFTER RO.O t>AYS 



Fig. 18. Balloon positions, on a polar equal area projection, at day 80 
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tSOPYCKEC SURFACE ao*0,32<» KG/K**3. StTUATIOK AFTER 85.0 DAYS 



Fig, 19. Balloon positions/ on a polar equal area projection/ at day 85 
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Fig. 20. Balloon posiHons, on a piar equal area projecHon, ai- day 90. 
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1$09YCNIC $URfACt SlTlIATinH AFTER 95.0 JIAYS 



Ficj. 21 . 


Balloon posifions, on a polar equal area projection/ at day 95.. 
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ISOPVCNIC SlDtFACE Rn>0.3?9 StTU*T|nN ATTER 100.0 DAYS 



Fig. 22. Balloon positions, on a polar equal area projection, at day 100 
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Fig. 23. Balloon posUionS/ on a polar equal area prciviofion^ at day 105. 
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ISOpVcHIC ft0<0.32« KG/M««3« StTU^HON AFTER 110. 0 Q&rS 



Fig. 24. Balloon poslHons, on a polar equal area projection, at day 1 10. 




Fig. 25. Balloon positions, on a polar equal area projection, at day 115. 
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ISOPYCNfC StIftFACe RO»0.‘»29 SIIUATIDN AFTER 120.0 OATS 



Fig. 26. Balloon posltons^ on a polar equal area prolection, at day 120 
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liOPYCNIC SUKfACE Rft-O.WS Kr,/M**3, SlPIATfON APTER I?^,0 &*VS 



Fig. 27. Balloon positions, on a polar equal area projection, at day 125. 
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ISUrrCKIC SU'^F^CF KG/M**3. SjrtllTIOK iFTt« D*YS 



Fig, 28. Balloon positions/ on a polar equal area projection.- at day 130. 
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By about 30 days the balloons approach a quasi-steady state, in which 
almost all of the B and C balloons, and some of the A balloons, are within a band 
centered near 45°S latitude. Within this general zonal distribution there are several 
clusters. These clusters change in shape and intensity as they drift eastward with 
time, circling the globe in the order of ten days. The more intense and larger dusters 
can be traced on the consecutive 5 day maps. 

From day 70 to day 80 a large cluster forms near the Ross Sea. it splits 
into two smaller clusters by day 90^ and disperses almost completely by day 110. 

From about day 85 to day 120 a single large cluster, near 45°S latitude, 
is very conspicuous in the balloon distribution. After day 120 if weakens by zona! 
elongation. 


he gradual mixing of the three kinds of balloons can be followed on the 
consecutive maps; but it proceeds at a slow rate. Throughout the 1 30 days only a 
small number of the B and C balloons enter the central (polar) region. The central 
region remains mainly populated by the A balloons. On the other hand, a large 
fraction of the A balloons form part of the middle latitude zone of high density of 
balloons. 


After day 25 (26 December, year 1), no balloons move northward across 
the 30°S latitude circle in the numerical simulation. After day 75 very few balloons 
move northward across 35°S latitude. 
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5. STATISTICAL DESCRIPTION OF THE BALLOON DISTRIBUTIONS 

Two types of statistical description were made of the space-time distri- 
butions of the balloons In this simulation of the 1970-71 Eoie Experiment. 

In the first description, the degree of clustering Is shown by taking 1000 
points, located at random on the spherical earth poleward of 30°S latitude, and • 
computing the overage of the distance from each of these points to the nearest 
balloon. It was shown by Mesinger (1965) that if the balloons are randomly distri- 
buted with respect to each other, this average nearest balloon distance has a minimum 
value when the balloons have a constant density in space (a constant number per unit 
are). This minimum value is 0.5 when the unit of length is chosen to make the average 
density of balloons equal to 1 per unit area. (In the present study, when all 500 
balloons are poleward of 30^ latitude this unit of length is 505 kilometers; and the 
average nearest balloon distance, for balloons that are randomly distributed with a 
constant number per unit area, is 252,5 kilometers.) If the balloons deviate from the 
desired state of random distribution with a constant number per unit area, then the 
mean nearest balloon distance will be larger than 0.5; and its deviation from 0.5 will 
indicate the degree of balloon clustering. 

Figure 29 shows the average nearest balloon distance as a function of time, 
as computed ct 2.5 day intervals; starting from day 17.5, when all 500 balloons are 
in the air. 

At day 17.5 the nearest balloon distance is 0.83 (approximately 420 km). 

It slowly increases vrith time to its maximum of 1 .30 (657 km) at day 105. 



NEAREST BALLOON DISTANCE 
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Fig. 29. Average distance to the nearest balloon from 1000 points located 
at random poleward of 30°S latitude (at intervals of 2.5 days). 

The unit of length is chosen to make the average density of balloons 
(the number per unit area)/ in the region poleward of 30 S , equal 
to 1 . (When 500 balloons are poleward of 30 5/ this unit of length 
is 505 kilometers.) 

if we compare the map for day 105 (figure 23) with, say, the map for 
day 20 (figure 6), we see that on day 105 there are several strong clusters of 
balloons and also large empty areas near the 30° latitude limit. These strong 
clusters accompanied by large empty areas are not pronounced on the map for 
day 20. Day 85 has the smallest average nearest balloon distance during the 
last part of the simulation. The map for day 85 (figure 19) shows a smaller number 
of clustei's than day 105 and smaller empty areas near the 30° latitude limit used 
in calculating this index. 
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In the second description of the balioon distributions, in figure 30, we 

4 

show the distribution of the balloons with latitude, as averaged for the last 50 days 
of the simulation (day 80 to day 130). The last 50 days were chosen because the 
time-variation, shown in figure 29, indicates that the balloon distribution may 
already be in a statistically steady state. Figure 30 shows the average of 21 
histograms, taken at 2.5 day intervals, in each of which the average density of 
balloons is computed for 5^ latitude intervals. The ordinate scale on the left shows 
the density of balloons (the number per unit area) when the unit length is chosen 
such that the average density over the entire globe is 1 . The ordinate scale on the 
right corresponds to the scale used in figure 29, where the unit length was chosen 
such that the average density over the polar cap south of 30 S latitude is equal to 
in which case the unit length Is 505 km. The abscissa scale in figure 30 is propor- 
tional to the sine of the latitude. 

Figure 30 shows that there is more than one balloon per (505 krn) in the 
latitude zone between 40*^S and 55°S;and more than two balloons per (505 km) 
in the zone 45°S to 50°S. There is a.secondary maximum in the polar region, south 
of 75°S latitude. From 40°S to 35°S there is an average density of about 0.5 
balloons per (505 km)^. There were only three occasions (In the 21 histograms) in 
which a single balloon was between 35 S and 30 S. As inciated earlier, in none 
of these 21 histograms was there a balioon equatorward of 30 S latitude. 

The distribution shown in figure 30 is based on a large number 
(500 X 21 = 10,500) of balioon locations. The two maxima of balloon density, 
and the absence of balloons equatorward of 30 S latitude indicate the upper 
tropospheric convergence and divergence of the three cell mean meridional 
circulation produced in the numerical simulation of the global circulation. 



DENSITY OF BALLOONS 
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■ LATITUDE, degrees S 


Fig. 30, The average density of balloons (the number per unit area) In each 
five degree latitude zone, averaged for the last 50 days of the 
numerical simulation. The abscissa scale is proportional to the 
sine of the latitude. The ordinate scale on the left gives the density 
of balloons when the unit of length Is chosen such that the average 
density over the entire globe is 1 . The ordinate scale on the right 
shows the number of balloons per (505 km)^. The dotted line repre~ 
sents the density when all the balloons are uniformly distributed 
poleward of 30°S latitude. 

. As indicated earlier, figures 1 through 28 show that very few of the B 
and C balloons (those released at Neuque'n and Mendoza) enter the polar region; 
whereas a large fraction of the A balloons (those released at Logo Faghano) move 
into middle latitudes. This.suggests that we should examine the space-time distri- 
bution of just the A balloons by themselves. 
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in figure 31 , fhe solid curve is from figure 29, and is based on a!i 
500 A, B and C balloons. The dashed curve is fhe same statistical representation 
based only upon the 167 A balloons (the Logo Fdgnano balloons). From about 
day 30 on, the average of the nearest balloon distance, for the A balloons alone, 
is closer to the value for a random balloon distribution (0.5), than is the curve 
for the balloons from all three launch sites. This, of- course, is only the statistical 
statement of what the maps showed - 



TIME, days 


Fig. 31 . Average distance to the nearest balloon from 1000 points located 
at random poleward of 30 S latitude (at intervals of 2.5 days). 

The unit of length is chosen to make the average density of balloons 
(the number per unit area), in the region poleward of 30 S, equal 
to 1 . The solid curve is the distribution shown in figure 29, based 
on balloons A, B and C. The dashed curve is the distribution 
obtained with only the 167 A balloons (those released at Logo 
Fagnano). 
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In figure 32, the solid curve is from figure 30, and is based on all 
500 A, B and C balloons. The dashed curve Is the same statistical representation 
for the 167 A balloons only. The more pronounced polar maximum of the A balloons 
alone, reflects the fact that relatively few B and C balloons entered this region. 
From 70°S to 45°S the distribution of the A balloons is more uniform with latitude 
than the distribution which included the B and C balloons. But between 40 S and 
35°S the density of the A balloons alone is poorer than the density for the A, B and 
C balloons. 


6. TEST COMPARISON WITH OBSERVED GHOST BALLOONS 


To evaluate the reliability of the results given above, we made a 
test comparison v/ith observed GHOST balloon trajectories. 

Solot (1968) showed the trajectories of 9 balloons at the density surface 
(for 200 mb). Of these, 7 balloons were In the air during the 100 day period 
1 January to 1 1 April 1967. All the GHOST balloon data for this 100 day period 
was used in our test comparison. Table I shows the time period and initial location 
of the GHOST balloons. 
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LATITUDE, degrees S 


Fig. 32. The average density of baHoons (the number per unit area) in each 
five degree latitude zone, averaged ibr the last 50 days of the 
numerical simulation experiment. The abscissa scale and the left 
side ordinate scale are the same as in figure 30. The solid curve 
is the distribution shown in figure 30, based on balloons A, B and 
C. The dashed curve is the distribution obtained with only the 167 A 
balloons (those released at Lago Fagnano). The dotted line represents 
the density of balloons when they are all uniformly distributed pole- 
ward of 30*^S latitude. 
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TABLE I 

GHOST BALLOONS (Data from Solot/ 1 968) 


Balloon Letter 

NCAR 

Identification 

Period 

Initial Location 

A 

49R 

1 Jan - 1 2 Jan 

31.55°S/ 124.71°W 

B 

34J 

1 Jan - 1 3 Jan 

46.43°S, 134.91°E 

C 

50D 

1 Jan - 1 Feb 

64.03°S-, 44.72°W 

D 

58G 

1 Jan - 25 Feb 

53.21°S, 62.20°E 

E 

60TA 

1 Jan - 1 1 April 

59.40°$/ 46.I2°W 

F 

70K 

19 Jan “ n April 

55.94°S, 160.41°W 

G 

79R 

17 Feb - 11 April 

45.97°S/ 148.77°W 


There were 5 GHOST feaUoons In the air on January Isf^^ 1967 , Their 0.0 GMT 
positions on 1 January are therefore taken as the initial locations for this study. 
Two balloons were launched latar^ and their first locations are given by Solot for 
19 January and 17 February respectively. Only balloons E, F and G remained in 
the air until the end of our 100 day test period. During the 100 day test period 
there were never less than three nor more than five GHOST balloons in the air at 
any one time. 

Taking the once per day (0.0 GMT) position of each of these GHOST 
balloonSf we obtain a total of 351 positions/ which are shown in figure 33. This 
gives an average of 3.5 balloon positions per day. 

At each intitial location and Initial time of the 7 GHOST bal loons, 
shown in table I, we introduced into the numerical simulation of the global atmos- 
pheric circulation/ at p 2 (for 200 mb), a simulated GHOST balloon, each of which 
was allowed to stay in the air for the same period as the corresponding true GHOST 
balloon. The once per day (O'.O GMT) positions of the simulafed balloons, which 
also total 351 positions, are shown in figure 34. 
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!t Is important to stress two factors v/hen comparing the true and the 
simulated GHOST balloon positions. The first is that the year of the true GHOST 
balloon trajectories is year 1 967 AD , whereas the year of the simulated GHOST 
balloon trajeciories in the year ^ in the numerical simulation of the atmosphere. 
Only the initial calendar day and the initial location on that day, were the sajjie 
for each balloon, A through G, 

The second important factor, in the test comparison, is that the number 
of balloons, never more than 5 atony one time. Is a very small sample. 

Comparing figures 33and 34weare impressed by their similarity. In 
both cases the balloons are uniformly distributed in the region poleward of 45*^S 
latitude, with a rapid decrease of balloon density equatorward of that latitude. 

In this time“space representation, a cluster is produced if a weak wind 
allows a single balloon to stagnate in a given region: But in the preceding space 
representations, figures 1 through 28, clustering was produced when different 
bal loons v/ere brought together by the convergent component of the wind field. 

Figure 35 shows histograms for the 351 locations shown in the maps of 
figures 33'and 34, We notice the absence of the middle latitude maximum that 
was seen in figure 30. By comparing the maps of days 30 to 80 (figures 8 to 18) 
with the maps of days 80 to 130 (figures 18 to 28), we see that the difference in 
' histograms cannot be accounted for by any difference of the circulation regime 
between the first half and the second half of the 1 00 day period used to obtain 
the lower histogram in figure 35. Therefore, the difference in histograms must be 
due to the insufficient size of the GHOST defined sample. The small sample also 
explains the large magnitude of the polar maximum in the simulated GHOST distri- 
bution. This maximum is produced by a single balloon O'calloon E) stagnating in 
the zone between 70° and 80°S latitude. 
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Fig. 33. Once a day posiHons, on an equal area proiectlon, of all the obser v ed 
GHOST balloons which were in the air during the first 100 days of 
1967. (See table I and text for identification of the balloons A ihrough 

G.) 
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Fig. 35. Observed and simuiated average densiiy of fhe GHOST balloons 

(ihe number per unif area) in each five degree iafltude zone, computed 
from the balloon positions shown in figures 33 ond 34. The ordinate 
scale is the same as on the left side of figure 3U. 
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On the other hand^ we see in figure 33 and in figure 35 that one of the 
true GHOST balloons was in the latitude zone 30° to 20°S on four days; but that 
this was not the case for any of the simulated GHOST balloons. Although year 
1 967 and year 2 constitute a very small sample of years, this discrepancy may 
indicate a true, though small, difference between the real atmosphere and the 
numerically simulated atmosphere (in the sense that the real atmosphere may have 
a somewhat smaller constraint on the migration of balloons into the upper troposphere 
tropics). 


On the whole, the overall similarity between the observed and simulated 
GHOST balloon distributions in this test comparison Is close enough to make us 
believe that the space-time distributions of the large number of balloons in the 
planned 1970-71 Eole Experiment will be approximately the same as the simulated 
distributions shown in sections 4 and 5 above. 


7. SUMMARY AND CONCLUSIONS 


In our numerical simulation of the Eole Experiment, nearly ail of the 
balloons end up in the region poleward of 35°S latitude, with a distribution in which 
the average distance from random points to the nearest balloon is about twice as 
large as it would be if the balloons were distributed at random within that region. 

In the numerical simulation, there is a maximum density of balloons In a 
zone near 45°S latitude, where the density is about twice the average value. 
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Within this zone, in turn, there are usually several clusters. These clusters change 
in shape and intensity as they move eastward around the globe, making one revolution 
in the order of 10 days'. Another, but less intense maximum density of balloons is 
found in the polar region, where clusters occasionally form. 

A test comparison, with the small sample of 7 observed GHOST balloons, 
indicates that in the real atmosphere of the Southern Hemisphere summer the balloons 
move closer to the equator than in the simulation experiment. But in all other aspects 
the test comparison indicates a high reliability of the numerical simulation. 

The fact that very few of the balloons launched at Mendoza and Neuquen 
enter the polar region, whereas the balloons launched at Logo Fagnano not only are 
well dispersed within the polar region but also have a fairly good latitudinal distri- 
bution in the middle latitudes, suggests that, with these particular launching sites, 
one would obtain a better {more homogeneous) space distribution of the balloons by 
launching more than one third of the total number of balloons at Logo Fagnano. 
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